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ABSTRACT 

We present Spitzer 3.6 and 4.5 fim photometry and positions for a sample 
of 1510 brown dwarf candidates identified by the WISE All-Sky Survey. Of 
these, 166 have been spectroscopically classified as objects with spectral types 
M(l), L(7), T(146), and Y(12); Sixteen other objects are non- (sub) stellar in 
nature. The remainder are most likely distant L and T dwarfs lacking spec- 
troscopic verification, other Y dwarf candidates still awaiting follow-up, and 
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assorted other objects whose Spitzer photometry reveals them to be back- 
ground sources. We present a catalog of Spitzer photometry for all astrophys- 
ical sources identified in these fields and use this catalog to identify 7 fainter 
(4.5 /iin ~ 17.0 mag) brown dwarf candidates, which are possibly wide- field 
companions to the original WISE sources. To test this hypothesis, we use a 
sample of 919 Spitzer observations around WISE-selected high-redshift hyper- 
luminous infrared galaxy (HyLIRG) candidates. For this control sample we 
find another 6 brown dwarf candidates, suggesting that the 7 companion can- 
didates are not physically associated. In fact, only one of these 7 Spitzer brown 
dwarf candidates has a photometric distance estimate consistent with being 
a companion to the WISE brown dwarf candidate. Other than this there is 
no evidence for any widely separated (> 20 AU) ultra-cool binaries. As an 
adjunct to this paper, we make available a source catalog of ~ 7.33 xlO 5 
objects detected in all of these Spitzer follow-up fields for use by the astro- 
nomical community. The complete catalog includes the Spitzer 3.6 and 4.5 
(Urn photometry, along with positionally matched B and R photometry from 
USNO-B; J, H, and K s photometry from 2MASS; and Wl, W2, W3, and 
W4 photometry from the WISE all-sky catalog. 

Subject headings: stars: brown dwarfs - stars: photometry - galaxies: pho- 
tometry - galaxies: high redshift 



1. Introduction 



In the historical nomenclature of stellar classification, astronomers have classified 
stars according to their spectra, finally settling on the follo wing spectral types, OBAFGKM, 
with O representing the hottest stars and M the coolest (IMorgan et al.lll943[ ). With the 
discovery of numerous (sub)stellar objects cooler than the coldest M dwarfs, the need for 
new spectral types became evident. To describe and characterize this new population 
two new spectral t ypes were introduced to the stellar taxonomy, L and T dwarfs (see e.g. 
Kirkpatrickl ( 120051 ) for a n in depth review), while a ne w type, "Y" , i s used for even colder 
objects, as suggested by iKirkpatrick et al.l (119991 ) and iKirkpatrickl (J2000). 



These low-mass, low-temperature astrophysical sources are considered to be the 
lowest-mass population of st ar-like objects. Those with masses be low ~ 0.075 M are 
referred to as "brown dwarfs" ( lKumarlll963l ; lHayashi fc Nakanolll963l ). and have core tem- 
peratures insufficient to sustain hydrogen fusion. Brown dwarfs are important objects 
for several reasons. They provide benchmarks with which to te st our theories o f gas 
chemistry and cloud formation in low-temperature atmospheres ( IBurgasserl 120081 ). As 
analogs of extra-solar planets brown dwarfs are easier to study because they are typi- 
cally isolated. They represent in situ time capsules of star formation, since their mass 



-3- 



is never recycled back into the interstellar medium, the refore preserving information on 
metallicity enrichment over the lifetime of the Galaxy (IBurgasserl 12008) . They provide 
important constraints on both t he functional form of the stellar mass function and the 



low-mass limit of star formatio n (IKirkpatrick et al.ll2012l ; iMetchev et al.ll2008t iReyle et al. 



2010i ; iBurningham et al.ll2010l ). Though predicted to exist in the 1960's, brown dwarfs re- 
mained elusive for de cades until surveys such as the Two Micr on All- Sky Survey (2MASS; 



Skrutskie et al.ll2006h . the Sloan Digital Sky Survey (SDSS; 



York et al. 



2000), and the 
19971 ). revealed 



Deep Near- Infrared Survey of the Southern Sky (DENIS; lEpchtein et al. 
them in large numbers. 

Recently researchers have shifted their focus toward the coldest brown dwarfs, which 
are one of the two primary science objectives of the Wide-field Infrared Survey Explorer 
(WISE; IWright et al.l l2010h . WISE launched on 2009 December 14 and began surveying 
the sky on 2010 January 7 at wavelengths of 3.4, 4.6, 12, and 22 /xm, hereafter referred to 
as Wl, W2, W3, and W4, respectively. It completed its first full pass of the sky on 2010 
July 17 and the WISE All-Sky Data Release was issued on 2012 March 140. The survey 
reaches 5a point source sensitivities in unco nfused regions to better than 0.08, 0.11, 1 
and 6 mJy, respectively ( IWright et al.ll2010l ). The Wl and W2 filters were specifically 
designed to probe the deep 3.3 jttm CH 4 absorption band in the spectra of brown dwarfs 
and the region relatively free of opacity at ~ 4.6 /xm, making their Wl — W2 colors very 
red and allowing cool brown dwarfs to be readily identifiable. 

Hundreds of potentially ultra-cool brown dwarf candidates have been identified using 
the W ISE All-Sky Data. A detailed description of this search is described in 



(2011) and Kirk patrick et al. f|2012h. where, alo n g with iMainzer et 



(120 llh . iBurgasser et al.l (1201 If ). ICushing et all (I201lh . iGelino et al 



(201 



Kirkpat rick et al 



UWright et al 



( 1201 lh . and Mace et 



al. (in prep.), follow-up and preliminary analyses have been presented. With the relatively 
shallow depth of WISE, a majority of these ultra-cool brown dwarf candidates either have 
very faint Wl detections or were simply not detected in Wl, resulting in Wl — W2 color 



limit s. The Spitzer Infrared Array Camera (IRAC) 3.6 and 4.5 /iin bands (IFazio et al. 



2004), sometimes referred to as chl and ch2 are very similar to the Wl and W2 filters in 
WISE. Thi s makes IRAC an ex cellent instrument for deeper follow-up observations (see 
Figure 2 in lMainzer et al.l 1201 if ). 



In this paper, we present the photometric properties from IRAC follow-up of a large 
sample (~1500) of WISE-selected brown dwarf candidates (see Section 4.1). Of these, 182 
have been spectroscopically classified, with classification s ranging from a single M dwarf 



to 12 examples of the latest an d coldest spectral type, Y (lKirkpatrickl|2008t ICushing et al. 



20 lit Kirkpatrick et al.ll2012l ). We also present new, ultra-cool brown dwarf candidates 
discovered using the Spitzer data alone in an effort to identify fainter, widely separated 
companions to the WISE sources. As a control sample we also identify field brown dwarf 
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-4- 



candidates from a Spitzer campaign to follow-up WISE-selected hyper-luminous infrared 
galaxy (HyLIRG) candidates (see Section 4.2), as this enables us to assess whether the 
Spitzer-selected brown dwarf companion candidates are more likely to be truly associated 
with the WISE-selected brown dwarf candidates. We release a catalog for 906 WISE- 
selected HyLIRG candidates. In addition to these two extremely rare population of as- 
trophysical sources, we release a photometric catalog of ~ 7.33 xlO 5 sources detected in 
the 1000's of Spitzer Astronomical Observation Requests (AOR's) in these follow-up cam- 
paigns. In §2 we briefly describe the observations. We discuss the photometry and source 
detection in §3. In §4 we describe the catalogs presented in this release. We discuss the 
new brown dwarf candidates discovered by this analysis in §5 and summarize the paper 
in §6. All magnitudes are given in the Vega system unless otherwise noted. 



Observations 



Warm Spitzer observations were carried out at 3.6 and 4.5 ixm under Spitzer programs 
70062, 70162, 80033, and 80109. Both the brown dwarf and HyLIRG candidate field 
observations used a 5-point dither pattern with 30 second exposures per pointing in each 
IRAC band. IRAC has a 5'x5' field of view with 1'.'2 pixels. For this analysis we utilized 
the post-BCD Spitzer pipeline images, which have been resampled onto 0.6" pixels. The 
AORs were executed between June 2010 and May 2012 and comprise a total of 1564 
brown dwarf and 919 extragalactic follow-up o bser yations. These observati ons have been 
described in detail in iKirkpatrick et al.l (120 llf ) and lEisenhardt et al.l (120121 ) . The images 
for this analysis were reduced with two different versions of the Spitzer reduction pipeline, 
18.18 and 19.1, with the only change between the two versions being an update to the 
masking process which includes knowledge of latents going back to previous AOR's. None 
of these changes affect our reduction and processing techniques. 



3. Detection and Photometry 



Source detection and photometry were carried out using SExtractor fjBertin &: Arnouts 



1996). We constructed a ch2 selected catalog by using the dual-image mode capabilities of 
SExtractor. In dual-image mode, sources are detected and their centroids and apertures 
are defined in one image and subsequently the photometry is m easured in anothe r imag e 



using those predetermined apertures and source centers. As in lEisenhardt et al.l ( 12010I ). 
the choice of using ch2 to construct this catalog was motivated by the fact that cool 
brown dwarfs are generally brighter in ch2 than chl. Photometry was measured in 6'.'0 
diameter apertures and corrected using the calibration values determined by the IRAC 
instrument team, -0.133 and -0.113 for chl and ch2, respectively. We use the IRAC zero 
points provided by the Spitzer Science Center (SSC), which are 18.789 and 18.316 Vega 
magnitudes for chl and ch2, respectively. SExtractor was configured to define a source as 
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a set of twenty or more connected pixels, each lying 1.0a above the background. We utilize 
the post-BCD coverage maps to provide an exposure flag (COVCH^) for all sources in 
the full catalog, with COVCH^ equal to the number of frames going into the post-BCD 
stack at that position. We also provide the SExtractor parameter CLASS_STAR to aid in 
separating point sources from extended sources. We compared the SExtractor aperture 
magnitudes for our brown dwarf candidates to an independent measurement of aperture 
photometry as presented in Kirkpatrick et al. (2011,2012). This independent photometry 
tool was written in IDL using both public scripts from IDL Astronomy User's Library 2 ] 
and proprietary code created specifically for this task. In short, this IDL code centroids 
on the known position of the brown dwarf candidate, obtains the aperture photometry of 
the source, and applies the appropriate aperture corrections as provided by the SSC. We 
find a mean offset Ach2 ~0.001 mag and a standard deviation sigma (Ach2) ~0.04 mag 
between the two codes, indicating that both are behaving reasonably. 



4. The Catalogs 



There are three separate catalogs presented in this paper. The first, described in 
§4.1 and Table 1, is a list of Spitzer and WISE photometry for 1510 WISE-selected brown 
dwarf candidates targeted in Spitzer programs 70062, 80109. The second, described in 
§4.2 and Table 2, is a list of Spitzer and WI SE photometry for 906 W ISE-selected HyLIRG 
candidates (hereafter, "WlW2-dropouts," lEisenhardt et al.l 120121 ) that were the targets 
for Spitzer program 70162 and 80033. The third, described in §4.3 and Table 3, is a list of 
all (COV > 3) sources found on the Spitzer IRAC images for the three programs above. 



4.1. The Brown Dwarf Catalog 



Brown dwarf candidat es were selected f r om th e WI SE source databas e s usin g the 
color criteria discussed in Kirkpatrick et al.l ( 120111 ) and iKirkpatrick et al.l (120121 ) . To 
summarize, our main selection used a color cut chosen to select cold brown dwarfs with 
types >T5. A color cut of W\ — W2 > 1.5 mag was used with WISE internal source 
lists from early processing runs and Wl — W2 > 2.0 mag was used for later processing, 
specifically using the combination of the WISE All-Sky Source Cat alog and WISE Al l-Sky 
Reject Table that are contained in the WISE All-Sky Data Release (jCutri et al.ll2012l ). We 
also performed another cut with a much more relaxed color criterion (W1 — W2 > 0.4 mag) 
to select bright (W2 S/N > 30), nearby brown dwarfs with types >L5. Other constraints, 
including the lack of a positional match in earlier-epoch all-sky surveys (2MASS, DSS1, 
DSS2), cuts on the W2 — W3 color, stipulations on the reduced x 2 value from Point 



2 http: //idlastro. gsfc.nasa.gov/ 
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Spread Function photometric fitting, etc. were also imposed, as discussed in Kirkpatrick 
et al. (2011) and Kirkpatrick et al. (2012). 

Most of the selected sources were scheduled for Spitzer chl and ch2 observations as 
part of programs 70062 and 80109. As of late-May 2012, a total of 1564 unique candidates 
had been observed in both IRAC bands. We were not able to identify a Spitzer counterpart 
for 54 of these candidates, which were therefore deemed to be spurious sources in WISE. 
Further inspection revealed that these spurious sources fell into two categories: (a) cosmic 
rays that bled through to the final images in the WISE coaddition process due to multiple 
cosmic ray hits and low frame coverage, and (b) WISE sources within nebulous patches (or 
nebulae themselves) rather than nearby brown dwarfs. For the remaining 1510 candidates, 
however, Spitzer counterparts were easily detected in both chl and ch2 bands. These 1510 
sources comprise the WISE-selected Brown Dwarf Photometric Catalog, as described in 
Table 1. 

The WISE and Spitzer photometry for each source in this catalog is graphically 
illustrated in Figure 1. This figure plots the Spitzer/IRAC chl — ch2 color as a func- 
tion of WISE Wl — W2 color. The 1510 sources from the catal og are shown by blac k 
dots. A total of 1 66 o f these have spectroscopic confirmation fromlMainzer et al.l (2011 ). 



Burg asser et al.l (120111 ) , Kirkpatrick et al.l (120111 ) , ICushing et al.l (120111 ) , Kirkpatrick et al. 



(120121 ) and Mace et al. (in prep), and are plotted with colored symbols, as indicated in 
the legend of the plot. The coldest brown dwarfs, the Y dwarfs, populate the reddest 
locus. The location of our latest type Y dwarf, WISE J182831. 08+265037. 7, which was 
targeted in a different Spitzer program (program 55 1; Mainzer, PI), is show n by the ma- 



genta octagon. As further illustrated in Figure 8 of Kirkpatrick et al.l (120121 ). the Spitzer 
chl —ch2 color for WISE 1828+2650, typed as >Y2, is not as red as objects typed slightly 
earlier, at Yl. 

Many excellent brown dwarf candidates in the catalog and on this diagram still lack 
spectroscopic data. Of the five objects in Figure 1 with colors of chl — ch2 > 2.8 mag, 
three are promising Y dwarf candidates - WISE J064723.23-623235.5 (chl - ch2 = 2.87, 
Wl - W2 > 3.77), WISE J082507.35+280548.5 (chl - ch2 = 2.96, W1-W2 > 3.75), 
and WISE J220905. 73+271 144.00 (chl - ch2 = 2.98, Wl - W2 > 3.68). All three are 
extremely faint in the J and H passbands, so we are awaiting spectroscopic confirmation 
using the WFC3 grisms onboard the Hubble Space Telescope. The other two sources, WISE 
035358.23+375458.5 (chl-ch2 = 3.00, W1-W2 = 3.52) and WISE J141127.45-612925.6 
(chl — c/i2=3.42, Wl — W2 > 3.35), are located within nebular patches on the Spitzer 
images and appear to be reddened or embedded objects rather than brown dwarfs. 

These last two sources highlight an important point: spectroscopic confirmation of 
sources is vital to establishing candidates as bona fide brown dwarfs. Sixteen of our 
targets, shown by the yellow dots in Figure 1, have been shown spectroscopically not to 
be brown dwarfs (Mace et al. in prep). These scatter widely over the diagram, proving 
that there is contamination at all colors. One interesting feature can be seen in the 
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lower left quadrant of the diagram where our L and early-T dwarf candidates lie. These 
objects have bright WISE and Spitzer magnitudes, so the colors of these objects are well 
measured. These high-quality measurements enable us to split objects with Wl — W2 < 
1.5 mag and chl — ch2 < 1.0 mag into two tracks. The lower track is replete with 
spectroscopically confirmed L and early-T dwarfs and primarily have 2MASS detections; 
the upper track has fewer confirming observations and generally no 2MASS detections, 
but objects with spectra here are typically not found to be brown dwarfs. Users of this 
catalog are cautioned that the sources here are candidates only, our selections were fairly 
liberal, so contamination by other types of sources is expected. 



4.2. The WlW2-Dropout Catalog 



Eisenhardt et al.l (120121 ) have identified a rare population of objects termed U W1W2- 
dropout" galaxies, because they are faint or undetected in Wl and W2 but well detected 
i n W3 or WA. Optical s pectroscopy shows that most of these objects have redshift z > 1.6 
(jEisenhardt et al.ll2012l ). which in combination with sub- millimeter foll ow-up detection s 
results in bolometric luminosities > 1O 13 L and in some cases > 1O 14 L (IWu et al.ll2012l ). 
This qualifies them as hyper-luminous infrared galaxies or HyLIRGs, the other primary 
science objective for WISE (the first being cool brown dwarfs, as noted in the intro- 
duction). The sub-millimeter data show their luminosities are dominated by dust with 
temperatures m ore than twice as high as other infrared luminous populations, which leads 
Wu et al.l (120 121 ) to refer to them as hot dust obscured galaxies or "hot DOGs" . A closely 
related population shows a high incidence (1/3) of very extended (> 30 kpc) Ly-a emis- 
sion or Ly-a blobs, which normally are found in less than 1% of galaxies, as described in 
Bridge et al.l (120121 ) . The extreme luminosities, hot dust, and prevalence of Ly-a blobs 
in this population suggests they may be a rare phase in the co-evolution of galaxies and 
their central supermassive black holes. 

Approximately 1,000 VFlVF2-dropouts have been identified in the WISE catalog, and 
because of their faintness in Wl and W2, 919 have been followed up in the Spitzer IRAC 
chl and ch2 bands. Of the 919 sources, 910 were selected in Spitzer cycle 7, 9 were added 
in cycle 8, and 13 failed to meet the SExtractor detection thresholds described in section 
3. 

The IRAC chl and ch2 bands are similar in wavelength to Wl and W2 but can 
easily reach much greater depths. These bands sample near-infrared wavelengths in the 
rest-frame at the typical z ~ 2 of V^H /F2-dropouts, providing a measure of the stellar 
population in the galaxies. Figure 1 of Eisenhardt et al.l (120121 ) shows the resulting dis- 
tribution of chl — ch2 colors converted to Wl — W2. Here we present (in Table 2) Spitzer 
and WISE photometry for the 906 sources detected in the Spitzer data by SExtractor. 



The bulk of the chl — ch2 colors lie well above Wl — W2 = 0.8, which indicates the 
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source has a high probability of being an Active Galactic Nucleus (AGN), as 



(120121 ) show. The optical spectra also often show AGN signatures (jWu et al. 



Stern et al. 



20121) , and 



it appears likely that the bulk of their extra ordinary luminosi t y is p owered by accretion 
onto supermassive black holes. As shown in lEisenhardt et al.l (120121 ) , the unusually high 
ratio of W3 and W4 emission relative to chl suggests a higher ratio of supermassive 
black hole mass to stellar mass than is found in local galaxies, implying the black holes 
grow before their stellar populations. This sequence is unexpected in the most frequently 



discussed scenarios of galaxy and AGN feedback (e.g. iHopkinsi 120121 ) 



4.3. The Complete Photometric Catalog 

We provide a catalog of the ~ 7.33 xlO 5 sources found in the Spitzer fields (COVCH1 
> 3 and COVCH2 > 3) whose properties are described in Table 3. We provide WISE and 
2MASS photometry from the WISE All-Sky source catalog for sources that fall within 
2" of a Spitzer source. In addition, we provide optical photometry for Spitzer sources 
that fall within 2" of a source in the USNO-B catalog. Table 4 summarizes the statistics 
that resulted from matching between these catalogs for two different samples. Sample A 
comprises sources with the deepest imaging and most reliable photometric measurements 
available (COVCH1 = 5 and COVCH2 = 5). Sample B represents sources with slightly 
shallower imaging (COVCH1 > 3 and COVCH2 > 3) and thus slightly lower photometric 
reliability, though still reasonably good enough for most scientific purposes. 

To characterize the photometric sample, we provide color-color and color-magnitude 
diagrams in Figure 2. The left plot of Figure 2 is similar to Figure 1, but includes all 
sources from Sample A having WISE measurements (a total of ~ 10 5 sources) rather 
than just the WISE brown dwarf candidates. We also show these colors for the 166 
spectroscopically confirmed brown dwarfs, plotted with colored symbols as indicated in 
the legend of the plot. Few field sources in the Spitzer images have colors as extreme as 
the latest T and Y dwarfs. 

In the center plot of Figure 2 we show chl — ch2 versus ch2 mag. The limiting 
magnitude is ch2 ~ 18.0 mag, which is ~ 2.5 mag fainter then the WISE-selected candi- 
dates and comparable to the depth of the IRAC Shallow Survey (Eisenhardt et al 2004), 
which had similar exposure time. We see that the majority of the WISE-selected brown 
dwarf candidates occupy a diffuse region in this color-magnitude domain, mainly with 
ch2 < 15.5 and chl - chl > 0.7. 

It is challenging to separate Galactic sources from extragalactic sources using the 
chl and ch2 filters in the IRAC camera alone. In the rightmost plot of Figure 2 we show 
chl — ch2 vs ch2 — W3. Here we can see the benefit of including photometric measurements 
at 12/im. The additional color information helps separate the Galactic and extragalactic 
red (chl — ch2 > 1) populations with cool brown dwarfs having ch2 — W3 < 3 while the 
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extragalactic population has ch2 — W3 > 3. A similar separation is shown in Figure 1 of 



Eisenhardt et all (1201(1 ). 



5. New Brown Dwarf Candidates Selected From Spitzer 

We present new brown dwarf candidates discovered from a blind search using the 
following search criteria chl — ch2 > 1.5, ch2 < 17.0 ma g and coverage equal to 5 in 



both IRAC channels (sample A from Table 4). As noted by lEisenhardt et al.l ( 120101 ) . this 
choice of color is expected for sources with spectral types later than T7 and the magnitude 
cut allows for only robustly measured sources to be inspected. This search resulted in 
666 candidates from the brown dwarf fields and 247 candidates from the extragalactic 
fields. We then used the Basic Calibrated Data (BCDs) to inspect the individual frames 
to remove any spurious sources or sources which appeared to be contaminated by cosmic 
rays. After also removing sources which were the WISE-selected targets, we were left with 
13 new ultra-cool brown dwarf candidates: seven in the WISE brown dwarf fields, and 
six in the WISE HyLIRG fields. Table 5 lists these 13 candidates, and Multi-wavelength 
finder charts (1' x 1') for them are presented in Figures 3, 4, and 5. The finder charts 
show that all of these candidates lack detections at shorter wavelengths and 12 /im, as 
expected for brown dwarf candidates. However, the Spitzer images are relatively deep 
compared to the observations at shorter (DSS,2MASS) and longer wavelengths (WISE). 
The difference in depth between these sets of observations makes it difficult to exclude 
certain types of extra-galactic objects such as AGN. 



For comparison lEisenhardt et al.l (120101 ) found 32 simi lar brown dwarf c andidates 



in the deeper Spitzer Deep Wide Field Survey (SDWFS) (lAshby et al.l 120091 ) covering 
10 square degrees, but selecting candidates to a fainter limit, ch2 < 18.5. Of these 32 
SDWFS brown dwarf candidates, 7 had ch2 < 17.0. After accounting for the coverage 
depth of 5 criterion, the 2483 IRAC pointings searched here cover a comparable area to 
SDWFS, suggesting a somewhat higher surface density of brown dwarfs candidates than 
in SDWFS. All of the SDWFS candidates with ch2 < 17.0 were among those rejected as 
probable AGN because they had ch2 — ch4 > 2. An analogous AGN rejection criterion 
using the ch2 — W3 color is not possible here because the W3 depth is not adequate, so 
it is plausible that many of the 13 brown dwarf candidates identified here are AGN. If so, 
it is possible that the somewhat higher surface density is attributable to AGN clustering. 

Determining whether any of the new candidates could potentially be widely separated 
companions can be achieved in two ways. First, we check the surface density of our brown 
dwarf candidates in the control sample versus that in the brown dwarf fields. We find a 
surface density of (1.5 ±0.6 deg -2 ) in the control sample and (1.0 ±0.4 deg -2 ) in the brown 
dwarf fields. The errors are calculated based on simple Poisson statistics. The fact that 
the control field has a larger surface density suggest that the candidates found are unlikely 
to be physically associated companions. Second, we can determine the spectrophotometric 
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distances to both the primary WISE-selected source and the n ew Spitzer-selected candi- 
date. We estimate the spectral type by using Figure 11 from iKirkpatrick et al.l (j201lh 
and the derived chl — ch2 color. We then estimate the photometric distances listed in 
Table 5 using the apparent magnitude and the absolute magnitude, where the latter is 
estimated using the abs olute magnitude versus spectral type relationship in Figure 13 of 
Kirkpatrick et al.l (l2012h . 



Most of the new candidates do not appear to be physically associated with the WISE- 
selected targets, but are more distant background sources. Assuming the WISE-selected 
candidates are indeed brown dwarfs, they reside at ~ 20 pc, while the majority of the 
new Spitzer-selected candidates are at ^ 40 pc. There is only one case for which both 
the primary and the candidate may be physically associated, J052439. 3+464631. 2, where 
both sources have estimated distances of ~ 20 pc. The estimated spectral types for these 
two sources are T6 and T9. The angular separation of 83.26" for this system corresponds 
to a physical separation of at least 1600 AU at 20 pc. 

The probability of discovering widely separated companions rests on the probability 
that such systems exist, which is itself a complicated function of mass, age, environment, 
and format ion mechanisms. An in depth review on multiplicity in very low mass systems 
is given in iBurgasser et al.l (120071 ) and references therein. These brown dwarf primaries 



have low mass and thus small gr avitational potential wells. This is well summarized in 
Figure 6 of lBurgasser et al.l (120071 ) . where they plot the separation of binaries as a function 
of total mass for a wide variety of known systems. Brown dwarf binaries tend to have 
small separations with higher binding energies. Systems with lower binding energies tend 
to be disrupted, primarily due to outside perturbers or gravitational instabilities, which 
translates to almost no old brown dwarf binaries having been found with separations 
greater than 20 AU. There are only a handful of low-mass binaries with separations 
greater than 20 AU, three are relatively young and thus hav e not had enough t ime for 
collisional disruption. Another system is DENIS 0551-4443AB ( IBilleres et al. l2005h . which 
is thought to be relatively old a nd is the kind of widely separated system for which we were 
searching. IBilleres et al.l ( 120051 ) suggest that such a system is fragile, and it would not have 
survived a close encounter with a third body. Furthermore, its existence demonstrates that 
some very low mass stars/brown dwarfs form without ejection from a multiple system, 
or any other strong dyna mical interaction . Oth er widely separated multiple systems 



have been discovered since Burgas ser et al.l (l2007f) . including the NLTT20346 + 2MASS 
J0850359+105716 system dFahertv et allboilh . Is it possible that J052439.3+464631.2 is 



also an old, widely separated binary like these? 

The nature of J052439. 3+464631. 2 can be probed further by considering the environ- 
ment, magnitudes, and colors of the individual objects in this system. Using the WISE 
All-Sky Image Archive we investigate a region ~ 5' on a side centered on the new can- 
didate and make the following observations. First, this system resides near the Galactic 
Plane which makes it more difficult to interpret given the high source density and complex 
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nebulosity. The pair resides at the edge of a 12^m bright nebula, thus suggesting that 
rather than a widely separated T6-T9 binary, these are more likely reddened or embedded 
objects and not brown dwarfs. Another possibility is that they could be physically asso- 
ciated T6 and T9 brown dwarfs, though not in a binary system, i.e. forming in the same 
stellar cocoon and around the same time. Spectroscopic observations of these candidates 
should be made to disentangle these scenarios. 

Other than J052439. 3+464631. 2, this search yielded no plausible widely separated 
companions, therefore suggesting that very low mass (>T7) brown dwarf binaries with 
separations of 20 to 2000 AU are extremely rare. 

6. Summary 

This paper summarizes initial Spitzer follow-up and analysis of some of the rarest 
astrophysical objects discovered to date. We present, for the first time, a very large sample 
of WISE-selected ultra-cool brown dwarf candidates. Of these, 184 have spectroscopic 
data with ~ 92% of them confirmed as bona fide ultra-cool brown dwarfs and of which 
~ 83% have spectral types later than T5. We also present 13 new, though more distant, 
ultra-cool (>T7) brown dwarf candidates discovered using the Spitzer data. Only one 
possible widely separated companion system was found, suggesting that widely separated, 
cold brown dwarf binaries are extremely rare. To date, the WISE team has only followed 
up a small percentage (~ 12%) of the total WISE-selected sample, providing follow-up 
opportunities by the astrophysical community on a host of interesting discoveries yet to 
be made. Given the wide range of (sub)stellar spectral types it probes, the sample will 
allow for a multitude of future statistical studies. In addition, this sample will provide 
excellent follow-up sources for the James Webb Space Telescope (JWST). 

There are several avenues of investigation yet to be fully explored which we summa- 
rize as follows. Detailed proper motion measurements can be made with the current data 
available (WISE + Spitzer), though long term monitoring will allow for full kinematics 
and robust distances to be determined. Several large parallax programs are currently 
underway monitoring a large number of the spectroscopically verified sources, though 
many remain to be monitored. Long term monitoring will also shed critical insight and 
constraints on variability in brown dwarfs, two examples being the study of atmospheric 
meteorology and possible detection of transits by planetary companions. Using the dif- 
ference between the WISE and Spitzer filter profiles will allow one to explore differences 
in the underlying spectra at these wavelengths. With no all-sky deep surveys planned at 
these wavelengths in the foreseeable future, this is the quintessential brown dwarf sample 
because it comprises the closest and brightest brown dwarfs possible. 

In addition to presenting these two very rare populations of astrophysical sources, we 
release a catalog of ~ 7.33 xlO 5 objects detected in the Spitzer observations, including 
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mult i- wavelength (2MASS, USNO-B, and WISE) photometry for a subsample. Providing 
multi-wavelength photometry for hundreds of thousands of astrophysical sources will allow 
for a number of scientific studies to be conducted, from fully characterizing the mid-IR 
properties of the coldest brown dwarfs to a more complete characterization of the most 
IR luminous galaxies in the universe. 

This publication makes use of data products from the Wide-Field Infrared Survey 
Explorer, which is a joint project of the University of California, Los Angeles, and the 
Jet Propulsion Laboratory/ California Institute of Technology, funded by the National 
Aeronautics and Space Administration (NASA). This publication also makes use of ob- 
servations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under a contract with NASA. This work 
is also based in part on observations made with the NASA/ESA Hubble Space Telescope, 
obtained at the Space Telescope Science Institute, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS 5-26555. Some 
of the spectroscopic classifications presented herein were obtained at the W. M. Keck 
Observatory, which is operated as a scientific partnership among the California Institute 
of Technology, the University of California and the National Aeronautics and Space Ad- 
ministration. This publication also makes use of data products from 2MASS. 2MASS is a 
joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center/ California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. This research has made use of the 
NASA/IPAC Infrared Science Archive (IRSA), which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with NASA. 
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Table 1: Brown Dwarf Photometric Catalog Parameters 



Parameter 


Example value 


Description 


DESIGNATION 


205628.91+145953.2 


WISE sexigesimal designation 


SPTYPE 


Y0 


Near-IR Spectral Type 


RA 


314.12046 


Right ascension in decimal degrees J2000 


DEC 


14.998111 


Declination in decimal degrees J2000 


CH1MAG 


15.873 


3.6 /im Vega magnitude 


CH1ERR 


0.017 


3.6 /im Vega magnitude error 


CH2MAG 


13.914 


4.5 /iin Vega magnitude 


CH2ERR 


0.003 


4.5 /im Vega magnitude error 


W1MPRO 


18.253 


WISE 3.4 /im profile-fit magnitude 


W1SIGMPRO 


0.000 


WISE 3.4 /im profile-fit magnitude error 


W1SNR 


-0.300 


WISE 3.4 /im signal-to-noise ratio 


W2MPRO 


13.928 


WISE 4.6 /im profile- fit magnitude 


W2SIGMPRO 


0.046 


WISE 4.6 /im profile- fit magnitude error 


W2SNR 


23.400 


WISE 4.6 /im signal-to-noise ratio 


W3MPRO 


12.003 


WISE 12.0 /xm profile-fit magnitude 


W3SIGMPRO 


0.270 


WISE 12.0 /im profile-fit magnitude error 


W3SNR 


4.000 


WISE 12.0 /im signal-to-noise ratio 


W4MPRO 


8.781 


WISE 22.0 /im profile-fit magnitude 


W4SIGMPRO 


NaN 


WISE 22.0 /im profile-fit magnitude error 


W4SNR 


0.800 


WISE 22.0 /im signal-to-noise ratio 



Note. — When W#SIGMPRO is either 0.000 or NaN it is the same as null in the WISE All-sky Source 
Catalog, see http://wise2.ipac.caltech.edu/docs/release/allsky 
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Table 2: WlW2-dropout Photometric Catalog Parameters 



Parameter 


Example value 


Description 


DESIGNATION 


000025 1+420708 5 


WISE sexiffesimal designation 


RA 


10475330 


Riffht ascension in decimal decrees ,12000 

-1 \_i± U CAj k7 V 1 -L k7 1 v_/ J. J_ 111 V i V. ' V . lllll 11 vi\j)_.l vV_/U T 7 V7 V7 V7 


DEC 


42.119005 


Declination in decimal decrees J2000 

J — S V. ' V.y 1111(1 V_J J. Vll 111 V I V. ' V.y 1111(11 vlv^t*! vVjU ' I \J \J \J 


CH1MAG 


16.759 


3 6 //m Vesa magnitude 

'..J t \J JUC/lll V v^U llliljL^lll U l^i. V^-L w 


CHI ERR 


0.039 


3.6 /im Vega magnitude error 


CH2MAG 


15.837 


4.5 fim Vega magnitude 


CH2ERR 


0.022 


4.5 /jm Vega magnitude error 


W1MPR0 


17.912 


WISE 3.4 /im profile-fit magnitude 


W1SIGMPRO 


0.241 


WISE 3.4 /im profile-fit magnitude error 


W1SNR 


4.500 


WISE 3.4 /im signal-to-noise ratio 


W2MPRO 


16.198 


WISE 4.6 /im profile-fit magnitude 


W2SIGMPRO 


0.176 


WISE 4.6 /im profile-fit magnitude error 


W2SNR 


6.200 


WISE 4.6 /im signal-to- noise ratio 


W3MPRO 


11.474 


WISE 12.0 /im profile-fit magnitude 


W3SIGMPRO 


0.114 


WISE 12.0 /im profile-fit magnitude error 


W3SNR 


9.500 


WISE 12.0 /im signal-to-noise ratio 


W4MPRO 


7.50 


WISE 22.0 /im profile-fit magnitude 


W4SIGMPRO 


0.092 


WISE 22.0 /im profile-fit magnitude error 


W4SNR 


11.800 


WISE 22.0 /im signal-to-noise ratio 



Note. — When W#SIGMPRO is either 0.000 or NaN it is the same as null in the WISE All-sky Source 
Catalog, see http://wise2.ipac.caltech.edu/docs/release/allsky 



-15- 



Table 3: Spitzer Photometric Catalog Parameters 



jr ai ameiei 


H V • 1 1 n tt 1 1 1 t ro 1 no 

r_jxampie vai ue 


1 locr'T'i'nfiA'ri 

LJxzbLL ipilUll 




-rA nsi Q/i 
r4Uoiy4oo 


Astronomical Observation Request 


O U rXV il/ 1 


DD 


Origin of Spitzcv follow-up (BE) or EXGAL) 






rugni ascension in uecimai uegiees ^jzuuuj 




-z ( .ooo4oy 


Declination in decimal degrees (J2000) 


PIT 1 A/TAP 


1/11 
14. loo 


3.6 /im Vega magnitude 


PU1 171315 


n nnQ 
U.UUo 


3.6 /im Vega magnitude error 


puoi\/r A P 

LvllZlVl/\(or 


14. 14 1 


4.5 /zm Vega magnitude 


PIIOITBB 
oiizii/rxrx 


n nnQ 

U.UUO 


4.5 /im Vega magnitude error 


PT A QQ QT^AB PTI1 
LLAoo.o lAlX-Olll 


n oQn 
u.yoU 


SExtractor star-galaxy separator measured 3.6 /im 


PT AQQ QTAR PIT 9 
UJ_/i\00_0 1 AIX_LvllZ 


u.you 


SExtractor star-galaxy separator measured at 4.5 /im 


PPVPII1 


K 



Total number of Spitzer frames at source position 3.6 /im 


PPA TP IT 9 


r 



Total number of Spitzeriidsn.es at source position 4.5 /im 


AAH MPBH 
VV HVlr txXJ 


1/11 QO 

14. iyz 


WISE 3.4 /xm profile-fit magnitude 


AAH QTPATPBPi 
VV loHorlvlr IXXJ 


n n^o 
u.uoz 


WISE 3.4 /xm profile-fit magnitude error 


AAH QWB 
VV ISlNrt 


qa nnn 
04.UUU 


WISE 3.4 /im signal-to-noise ratio 


AAf9T\TT3T?P 
VV Zlvlr sxKJ 


1 /i onQ 
14.ZUO 


WISE 4.6 /xm profile-fit magnitude 


AAf9QTPT\TPI30 
VV ZolLrlvlr sxKJ 


n n^R 

U.UOO 


WISE 4.6 /im profile-fit magnitude error 


AAf9QlM"R 


1 o /inn 
iy.4UU 


WISE 4.6 /im signal-to- noise ratio 


VV olvlr rUJ 


1 O fifil 
IZ.OOl 


WISE 12.0 /im profile-fit magnitude 


AAAQQTPA/TPT3P 


1ST a 1ST 


WISE 12.0 /im profile-fit magnitude error 


VV oolN It 


n ^nnn 
-u.ouuu 


WISE 12.0 /im signal-to-noise ratio 


VV41Vlr sxKJ 


y.uoou 


WISE 22.0 /im profile-fit magnitude 


W4SIGMPR0 


NaN 


WISE 22.0 /im profile-fit magnitude error 


W4SNR 


0.100 


WISE 22.0 /im signal-to-noise ratio 


JMAG 


15.304 


2 MASS J magnitude 


JERR 


0.055 


2MASS J magnitude error 


HMAG 


14.635 


2MASS H magnitude 


HERR 


0.063 


2MASS H magnitude error 


KMAG 


14.471 


2MASS K s magnitude 


KERR 


0.083 


2MASS K s magnitude error 


BMAG 


19.690 


USNO-B B2 magnitude 


RMAG 


18.370 


USNO-B R2 magnitude 



Note. — When W#SIGMPRO is either 0.000 or NaN it is the same as null in the WISE All-sky Source 
Catalog, see http://wise2.ipac.caltech.edu/docs/release/allsky 
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Table 4: Photometric Catalog Statistics 



Sample Coverage #w/Spitzer # w/WISE # w/2MASS # w/USNO-B 
A ^5 3.41xl0 5 l.OOxlO 5 3.90 xlO 4 7.75 xlO 4 

B >3 7.33xl0 5 2.02xl0 5 8.22 xlO 4 1.63 xlO 5 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I l.l I I I I I I I I I 



3 - 



a 2 

U 

i 

S3 

u 



1 - 



• M(l) 
*L(7) 

A < T5 (20) 
A > T5 (126) 
▲ Y (12) 
nonBD (16) 

# Mainzer 2011 



A 



- 




I I I I I I I I I 



2 3 
W1-W2 



Fig. 1. — Spitzer /IRAC chl — ch2 as a function of Wl — W2 color. Black points are 
observations of 1510 WISE-selected brown dwarfs detected by Spitzer . Those with spec- 
troscopic classifications are overplotted in color, described in the legend. Open symbols 
denote color limits in Wl — W2. We present the median uncertainties in each full mag- 
nitude bin of color for both axes and these are shown as the tick marks along the lower 
and left axes. The dotted line represents the 1-to-l relationship between the axes. 



Table 5: New Brown Dwarf Candidates 



Origin Candidate Candidate Primary Type a Type b d° est d d est Separation 6 



RA(J2000) DEC(J2000) chl chl - ch2 ch2 chl - ch2 pc pc 



BD 


02:39:59 


.8 


+70:40:20.4 


16 


.42+0. 


.04 


1 


.91+0. 


16 


14. 


.14+0. 


.01 





.78+0. 


.01 


T8 


T5 


40.6 


21.8 


36.27" 


EXGAL 


03:02:21 


.1 


02:16:09.8 


16 


,80+0. 


.05 


1 


.54+0. 


18 
























EXGAL 


04:12:09 


.3 


69:45:24.1 


16 


.94+0. 


06 


1 


.65+0. 


22 
























BD 


05:24:39. 


.3 


+46:46:31.2 


15. 


.57+0. 


02 


2 


.29+0. 


.11 


14. 


.52+0. 


.01 


1 


.03+0. 


01 


T9 


T6 


20.5 


21.5 


83.26" 


EXGAL 


08:21:51 


.3 


23:45:17.1 


16 


.63+0. 


.04 


2 


.45+0. 


.32 
























BD 


16:29:50. 


.0 


+25:45:25.9 


16 


.45+0. 


03 


1 


.68+0. 


.12 


15. 


.44+0. 


.01 


1 


.06+0. 


03 


T7.5 


T6 


44.5 


32.8 


139.13" 


BD 


16:53:42. 


.8 


-16:00:10.0 


16 


,91+0. 


.08 


1 


.69+0. 


33 


14. 


.61+0. 


,01 


1 


.00+0. 


.02 


T7.5 


T5.5 


54.92 


24.1 


44.98" 


EXGAL 


20:38:38 


.1 


-17:49:49.1 


15 


.32+0. 


01 


2 


.61+0. 


.12 






















h 


EXGAL 


21:41:45. 


.5 


01:20:14.5 


15. 


,95+0. 


.02 


2 


.27+0. 


19 
























EXGAL 


03:14:52. 


.0 


-53:52:41.9 


16 


,81+0. 


.04 


2 


.54+0. 


.40 
























BD 


03:04:26 


.6 


+43:31:10.9 


16 


.77+0. 


.05 


1 


.93+0. 


22 


14. 


.90+0 


.01 


1 


.22+0. 


02 


T8 


T7 


47.6 


23.0 


21.38" 


BD 


18:22:38. 


.6 


-67:50:03.9 


16 


.66+0. 


.04 


2 


.06+0. 


.22 


15. 


.12+0. 


.01 


1 


.09+0. 


02 


T9 


T6.5 


33.9 


26.7 


96.27" 


BD^ 


13:59:42 


.7 


-23:32:35.6 


16 


.56+0. 


.04 


1 


.61+0. 


.16 














T7 











a Estimated brown dwarf type for new candidate 

b Estimated brown dwarf type for primary 

c Estimated distance to the new candidate, in parsecs 

d Estimated distance to the primary, in parsecs 

e Distance between new candidate and the primary, in arc seconds 

/ primary Spitzer target for this observation was not detected in the imaging, thus no information is given for the primary target 
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W1-W2 CH1-CH2 CH2-W3 



Fig. 2. — Left: Spitzer/TRAC chl-ch2 as a function of WISE Wl - W2 color. Contours 
represent all sources with photometric measurements from these two instruments (~ 10 5 
sources). Objects with spectroscopic classification are shown by the colored points. Open 
symbols denote color limits in Wl — W2. Center: Spitzer /IRAC chl — ch2 as a function 
of ch2 magnitude. Contours show ~ 3.4 x 10 5 sources (sample A in Table 4). Objects with 
spectroscopic classification are shown by the colored points. We also show the 13 newly 
discovered candidates as yellow circles. Right: Spitzer /IRAC chl — ch2 as a function of 
ch2 — W3 color. Contours are similar to the left plot. Brown dwarfs with spectroscopic 
classification greater than T5 are shown in orange (T) and blue (Y). We also present the 
WISE-selected HyLIRG candidates in red. Open symbols denote color upper limits. The 
contours represent the 40, 60, 80, 90, 95, 97, and 99 percentiles of the distribution. Tick 
marks along the axes of the left panel are computed as discussed in Figure 1. 
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Fig. 3. — Multi- wavelength finder charts for the new ultra-cool brown dwarf candidates. 
We present lxl arcminute cutouts from various all-sky surveys alongside the Spitzer 
discovery images. From left to right, these cutouts are from the Digitized Sky Survey I 
band, 2MASS J and H, Spitzer chl and ch2 and WISE W3. The images are oriented 
with north up and east to the left. 
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Fig. 4. — Figure 3 (continued) 
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Fig. 5. — Figure 3 (continued) 
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